A scaling theory of quantum breakdown in solids 
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We propose a new scaling theory for general quantum breakdown phenomena. We show, taking 
Landau-Zener type breakdown as a particular example, that the breakdown phenomena can be 
viewed as a quantum phase transition for which the scaling theory is developed. The application of 
this new scaling theory to Zener type breakdown in Anderson insulators, and quantum quenching 
has been discussed. 
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The behavior of solid state electronic systems in non- 
equilibrium conditions, particularly the insulator- metal 
transitions leading to dielectric breakdown in constant 
electric field just beyond a particular threshold, has been 
a matter of considerable investigation. This started with 
Zener P, Q , considering the charge-excited carrier gen- 
eration process as tunneling through the band gap A, 
from ground state to excited states in one-dimensional 
band insulator, using the time-independent as well as 
time-dependent gauge 0, [1] ■ Recent years have noticed 
increasing interests in the nonlinear phenomena of quan- 
tum breakdown in many-body systems put in strong dc 
electric fields, in intense laser fields, etc. The interest also 
comes because of the fact that the high sensitivity (dras- 
tic changes) of electronic states to external field near the 
phase boundary finds its application in dielectric break- 
down based current-switching devices However, more 
interesting is its relevance to non-equilibrium statistical 
physics where one can observe the behavior of uncorre- 
latcd as well as strongly correlated electronic system near 
phase boundary under non-equilibrium conditions and in 
the study of quantum defect production in rapid quench 
of quantum fluctuations (see e.g., (H])- 

The tunneling rate P of charge excited carriers through 
a gap A in dc electric field E (for the case of band in- 
sulator) has been shown to be expressed in closed form, 
which is proportional to exp [—{£th/Ee)]; with Sth ^ A^ 
and e denoting the carrier charge. Recently, it has been 
shown that such Landau-Zener transition estimate can 
be straightforwardly extended to Mott insulators as well 
with the same exponent value 0, 0] , though the carriers 
are different. In band insulator holes and electrons un- 
dergo tunneling through band gap, whereas in Mott insu- 
lator doublons and holes carry out the same role through 
Mott gap. 

Within the present framework, it seems to be virtually 
impossible to obtain such an exact expression for the tun- 
neling rate of the carriers in a disordered quantum sys- 
tems like the Anderson insulators Q. Our purpose here 
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is to develop a simple scaling theory for such quantum 
breakdown phenomenon and explore the Landau-Zener 
type 3 breakdown in, e.g., Anderson insulators (say, in 
three dimension). It may be mentioned that the growth 
of fluctuation correlations near various well-studied clas- 
sical breakdown points have already helped establishing a 
phase transition picture, leading to well established scal- 
ing theories for phenomena, like fracture etc. (see e.g., 
0). Also, precise critical behavior of breaking in, e.g., 
the fiber bundle like models of fracture are now quite 
established [lo| . 

From elementary quantum mechanics (see e.g., ^3]), 
one gets the tunneling rate through an energy barrier of 
(large) width w to be proportional to exp(— ckw), where 
K is the damping factor or 'imaginary momentum' vec- 
tor determined by the height of the barrier and c is some 
dimensionless constant. Considering the quantum break- 
down phenomena as quantum phase transitions (see e.g., 
[l3|). we assume that near the breakdown point, the 
macroscopic correlation length ^ will provide the only 
scale governing the behavior of physical quantities of the 
system. We assume this to be true for the tunneling rate 
P as well. Hence, the width w is scaled by ^ and 'mo- 
mentum' is considered to be proportional to inverse of ^. 
Hence, the dimensionless tunneling rate expression takes 
the form 




Here, a denotes lattice constant, the only other (micro- 
scopic) length scale in the system. The above gives our 
proposed scaling form for tunneling probability across a 
barrier of width w near a quantum phase transition char- 
acterised by the correlation length ^ determined by the 
barrier height. 

We now proceed with the application of the above 
scaling ansatz. In presence of a constant dc field, the 
bands effectively get tilted in field direction, which causes 
each energy band degenerate with the others. Therefore, 
electron can pass from lower energy band to upper one 
traversing an effective distance w in the direction of field 
and we envisage the process involved as tunneling. The 
effective width w can be estimated assuming that the en- 
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ergy Eew acquired from the field E over the width w 
will compare with the gap A, giving w = A/Ee. We 
assume, the correlation length diverges as A^'' at criti- 
cality (A = 0) with exponent v [l^l, and additionally we 
assume that even up to the breakdown point the exter- 
nal field does not affect the correlation length exponent. 
Therefore the tunneling rate expression becomes 
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For the band and Mott insulators, putting f as 1/2 
(see e.g., [11]), we get P ~ exp{-£th/Ee),£th ~ A^, the 
same expression as was obtained earlier 0, 0, @, 01 • ^'^^ 
Mott insulators, the recent study by Oka et al. 0, 0] 
supports the above expression as well, where the gap is 
determined by the electronic correlation. 

Since the system response near breakdown is governed 
essentially by the correlation length exponent, one can 
make some estimate for the electrical response in quan- 
tum systems with disorder, equipped with the knowledge 
of its correlation length around criticality in such dis- 
ordered systems. Studies on Anderson transition show 
that the electron as a quantum particle can not diffuse 
through the geometrically percolating path due to the 
coherent back scattering (of the wave function) from the 
random geometry of the clusters in dimensions less than 
three Q . Since all the states on any such percolating lat- 
tice gets localized (exponentially), electrons do not dif- 
fuse through the disordered (classically percolating) lat- 
tice. In three dimension (for concentrations above the 
quantum percolation threshold) , if the Fermi level e f lies 
below the mobility edge Ec (f^ for three dimension), 
the system remains to be quantum mechanically non- 
percolating with the exponentially localized electronic 
states, while for e/ > ec the states are extended. Hence 
if £/ < the system is an insulator and for > Ec it 
is a conductor. The scaling property in this conducting 
region is well established; in particular the localization 
or correlation length C k/ ~ ^c\~'^ and precise esti- 
mates (both theoretical and experimental) of Ec and v 
are now available (see e.g. [1, 14|). For e/ < one 
can think of an insulator-metal transition in such a non- 
interacting random system, analogous to quantum break- 
down in band or Mott insulators, induced by (strong) 
dc electric fields. The critical behavior can be easily 
predicted utilizing the precise knowledge of localization 
length exponent v. With ~ 1 0] for Anderson insula- 
tors in three dimension, the tunneling probability across 
the mobility gap A = |ec — e/| near breakdown point, will 
be given by P ~ exp{—£th/ Ee), where £th ~ A'''; 7 ~ 3. 
This indicates that In P would scale as A^^ for Zener type 
breakdown in Anderson insulators, instead of as A^ for 
similar breakdown in band insulators. 

In the context of adiabatic quantum computations, 
where one exploits the tunneling probability through the 
energy or cost barriers for local minima to reach a global 



minimum (may be degenerate, as in spin glass), one needs 
to estimate the density of defects remaining over the 
ground state as one quenches from a highly excited (para) 
state [ij. For example, in a transverse Ising glass model 
represented by the hamiltonian Q 
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with random exchanges Jy between the Pauli spins <?, 
one can vary the quantum fluctuation slowly but linearly 
(r(t) = 1 — (t/r); < t < r) to arrive at a ground state 
of the classical spin glass (for F = at t = r) from a 
disordered or a para phase (for F = 1 at < = 0). The 
amount of defect over the true ground state (giving the 
solution of a computationally hard problem) depends on 
the quenching time r and the advantage of such annealing 
or quenching (tunneling) through macroscopically high 
or 0{N) barriers (instead of thermal hopping over the 
barriers as in classical simulated annealing [15[) might 
be lost if defect concentration remains high. Identifying 
the force Ee in Q as the rate R of change in the energy of 
the system, one can easily estimate the density of defects 
by calculating the tunneling probability as 
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where i? ~ ^ < H{t) >. Assuming that the residual en- 
ergy Eres (over the ground state energy) of the quenched 
state to be proportional to JA Pp{A) dA, and averag- 
ing over the Gaussian distribution p{A) ~ exp(— A^)[16[ 
of the local gap parameter A or of the local fields 
{hioc, A = \//j?oc + r^)[l|, one gets the decay rate of 
average residual energy as Eres ~ r~3 (for large t val- 
ues) in a long-range Gaussian spin glass model (having 

= 1/4 [l3|)- This, in fact, compares quite well with the 
numerical estimate of Eres for fairly larger Gaussian spin 
glass systems studied in [l8| . where Eres decays much 
slower than (observed for smaller systems). 

In brief, we propose here a generalized scaling form 
^ for the single particle tunneling probability across a 
barrier, where we subsequently replace the single particle 
correlation by the many body one appropriate for a quan- 
tum many body phase transition, obtaining thereby the 
generalized Landau-Zener breakdown probability. This 
scaling form has then been utilized for a Zener type 
breakdown in Anderson insulators in three dimension. 
Finally, we apply it to estimate the quenching rate de- 
pendence of the residual energy of a long range quantum 
Ising spin glass when its transverse field is swept through 
the critical point. Comparision with numerical results 
seems encouraging. 

We are thankful to Arnab Das, Jon Eriksen, Pradcep 
K. Mohanty and Purusattam Ray for important discus- 
sions. 
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